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For the 3-qubit UPB state, i.e^ the bound entangled state constructed from an Unextendable 
Product Basis of Bennett et al. , we provide a set of violations of Local Hidden Variable (LHV) 
models based on the particular type of reflection symmetry encoded in this state. The explicit 
nonlocal unitary operation needed to prepare the state from its reflected separable mixture of pure 
states is given, as well as a nonlocal one-parameter orbit of states with Positive Partial Transpositions 
(PPT) which swaps the entanglement between a state and its reflection twice during a period. 
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A bound entangled state is a nonseparable quantum 
state for which no distillation process is possible [2j. 
While some multiparty bound entangled states are eas- 
ier to detect because they show some form of bipartite 
entanglement through some of the cuts [3J, |^ , the bound 
entangled states built from an Unextendable Product Ba- 
sis (UPB states) of 0,0 are probably the most mysteri- 
ous one to date. For them, in fact, all partial transposi- 
tions are positive (PPT), no violation of Bell inequality is 
known and no criterion exist for detecting their entangle- 
ment. They are defined as the "complement" of a mix- 
ture (with all equal weights) of reciprocally orthogonal 
product states (such that no other product state orthog- 
onal to all members exists) in the given Hilbert space. In 
other words, since any separable (pure or mixed) state 
is a convex combination of pure product states [6j , if we 
take an orthogonal basis of such states which is unex- 
tandable, i.e., such that there does not exist any other 
product state orthogonal to the basis, then states which 
belong to the complement of this incomplete space in the 
Hilbert space of the quantum system are guaranteed to 
be nonseparable. If the mapping of the given mixture 
is done as in [l| mixing the UPB basis (with a minus 
sign in front) with the random state, then the property 
of PPT is also guaranteed. This mapping is nonunitary 
because of the minus sign just mentioned and it is shown 
in to be a particular case of a class of mirror-like 
symmetries whose origin and behavior are easily visible 
once we use the real tensorial parametrization proposed 
in 0. Such a reflection symmetry is nonequivalent to 
any known quantum symmetry and makes sense only for 
mixed states of multiqubits as the UPB states discussed 
here. It sheds considerable light into the structure of 
the bound entanglement of these states, as for instance 
it allows to identify violations of LHV models involving 
3 commuting observables, one violation for each term of 
the separable mixture that is reflected. We will see it 
for the 3-qubit UPB state (which is the quantum state 
considered in the entire paper). 

Two other issues are discussed after that: first the 



preparation of a UPB state, starting from the cor- 
responding separable mixture and using nonlocal uni- 
tary operations; second the existence of nonlocal orbits 
which arc PPT. We have found a single nontrivial one- 
parameter such orbit whose peculiarity is that the bound 
entanglement swaps between a state and its reflection 
twice during the orbit's period. 

Consider the 3-qubit density operator corresponding 

to the UPB state p UPB = { (l 8 - £}=i \ipj)(i>j\) , with 

IV*,-) = |01+),|1 + 0),| + Oil | ) (where |±) = 

-j= (|0) ± |1))) introduced in In terms of the tensor 

of coherences ^jj, if x — Pvpb is: 



3 = TT2 

{031, 033, 103, 111, 133, 303, 310, 313, 330, 331} 



PUPB 



Q 



{011,013,101,110,130,301} _ 



<? kl = 



otherwise. 



(1) 

The "complement" mentioned above is identified by the 
following operation: 
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PUPB = -Is — Pscp, 



(2) 



where p scp = 3X^=1 \ipj)(^j\- If we call p SC p 3 = iViXV'jl, 
in terms of the local Bloch vectors, p scp — j Pscp,, 
has components (g° — g 1 = g 3 = -^): 

Psc Pl = {g°Xo + g 3 X 3 ) <S> (g°Xo - g 3 X 3 ) <8> (g°\o + g 1 \i) , 
Psc P2 = (g°X - g 3 X 3 ) <8> (f?°A + g 1 Xi) ® (g°X + g 3 X 3 ) , 
Pso P3 = (g°X + ^Ai) <g> (g°X Q + g 3 X 3 ) ® (g°X - g 3 X 3 ) , 
Psc P4 = (£>°Ao - g 1 X 1 ) ® (g°X Q - g 1 Xi) ® (g°X - g 1 X 1 ) . 

(3) 

In pfl, it is shown that for the tensor of coherences the 
operation J2J is a reflection symmetry applied to the 
joint density of the 3 qubits, which results in a change 
of sign to the homogeneous part of the tensor, i.e., to 



2 



all the expectation values of 1-, 2- and 3-qubit observ- 
ables but not to the trivial one g 000 = tr (p SO pAooo) that 
"takes care" of the trace: if p sep = £> 000 A oo + £> then 
Pupb = S 6 4(p scp ) = £> 000 A 000 - £, in the notation of 0. 
The map 5*43 is a nonunitary linear map that preserves 
trace and Hermiticity. It is not well-defined in the set of 
admissible densities of 3 qubits, call it T>, although it is in 
the smaller subset C = {p £ V s.t. eig(p) £ [0, |]} where 
it is even invariant: p £ C implies p — Sq4 (p) £ C. Ob- 
viously p SO p £ C and pupb £ C. Because of the all equal 
weights, the mapping @ is also isospectral: eig(p scp ) = 
eig(puPB) = {0, j} of multiplicity 4. The reflection oper- 
ation (restricted to C) is obviously more general than the 
UPB construction of and needs not yield a PPT den- 
sity. For that the "unextendability" of UPB is required. 
No one of the 4 components p sePj taken alone (each is ob- 
viously a separable density) is a density when reflected. 
For example eig (S 6 4 (PsepJ) = {-§,3} of multiplicity 
resp. 1 and 7. Hence, although a posteriori (i.e., after 

the cancellations), it holds that C% J2^=i Pscp^j = 

\ Y^=i ^64 (psep^j "formally" the right hand side con- 
tains something else than just density operators. The 
convex statistical weights are fundamental, thus confirm- 
ing that operations like (0 are intrinsically defined only 
for mixed states. For the same reason the use and the 
meaning of the "complement" operation of must be 
handled with care. 

Look at the tensor of coherences of the 4 components 
of p sop as obtained from (J3J). It is straightforward to 
check that in the sum all and only the 1-qubit coher- 
ences are canceled. Such cancellations imply that all 
1-qubit reduced densities are random. A first conse- 
quence is that the same mapping p sep — > pupb can 
be obtained by means of two-qubit partial reflections: 
(<Si6 ® I4) (psep) = Pupb, and similarly for the other 
two. A second, more important, consequence is that 1- 
spin generators (— iadA j00 etc.) have no effect on pupb- 
Said otherwise, pupb is invariant to LOCC, even an in- 
finite amount of them, even stochastic, since the ran- 
dom reduced density is a fixed point also in the class of 
local filtering operations. A third consequence is that 
the 3-qubit bound entanglement cannot be modified by 
symmetric extensions of tensorial type 0, neither by 
convex combinations of extensions. From the expres- 
sion pupb = £ 000 Aooo — f 5 it is fairly easy to conclude 
on the character of any tensor product pupb ® Pa of 
Pupb with an ancilla p a . In fact, if for example the an- 
cilla is a qubit p a = g°\o + g{\j, then pupb <8> p a = 
e 0000 A ooo + ^"VaAoooj - £ <8> Q°Xo - £ <8> QiXj. Due to 
the affine structure, there is no way to modify the term 
— £ <8> g°Xo which is the "carrier" of the 3-qubit entangle- 
ment by means of LOCC (or local filtering). The result 
is the same for any ancilla, even entangled (even if we 
replace the 4-coherences term £ £g> g^Xj with something 
nontensorial) . 



Violations of LHV models. The UPB state is known 
to be bound entangled but no Bell inequality is avail- 
able for its detection. In we pointed out a class of 
multilinear algebraic inequalities which follow from the 
reflection symmetry and which lead to a simple form of 
contradiction of LHV models. Briefly (all details are in 
) such a contradiction has to do with the sign of the 
2-coherences of the reflection of a separable density in C. 
Unlike 1- and 3-cohercnccs, such sign pattern cannot be 
reproduced by simply changing sign in the Bloch vectors 
of the original convex sum because of a parity condition. 
If the original density has high rank and is close to the 
random state, then it is possible to compensate for the 
sign mismatch by adding a few more suitable terms in a 
mixture, but if the state is at the border of C and has 
rank 4, like pupb, the existence of a convex combination 
representing all 3-coherences of pupb and at the same 
time matching the "odd" sign pattern of its 2-coherences 
becomes unthinkable. This to justify why the inequal- 
ities reported below are not fully fledged entanglement 
criteria but rather contradictions, indicating a violation 
of a LHV model. Each of the 4 inequalities which can 
be obtained from (7b) of is corresponding to one of 
the 4 components of p sop in J3J, and each unveils a sign 
pattern which is not compatible with a LHV model: 

1) tr (pupbA 03 i) tr (pupbA 30 i) tr (pupbA 330 ) < 0, 

2) tr(pupBAoi 3 )tr(pupBAio3)tr(pupBAi3o) < 0, 

3) tr (pupbA 03 3) tr (pupbAi 3) tr (pupbAi 30 ) < 0, 

4) tr (pupbAoh) tr (pupbAioi) tr (pupbAho) < 0. 

Consider for example the case 1). Assuming the ex- 
istence of a LHV model, the measure "+" along Ao 3 i 
(we are only interested in its sign) must be compati- 
ble with the local measures fj(Xk) carried out on the 
j-th qubit along the A& axis. Two cases are possible: 
i) / 2 (A 3 ) - "+", h(Xi) = "+"; «) / 2 (A 3 ) = 
/s(Ai) = "— ". Similarly, along A301 we measure "— ", 
hence i) => /i(A 3 ) = it) / X (A 3 ) = "+". Since 

A 3 3o has measure "+", both the LHV models i) and ii) 
are contradictory. The inequalities are just a com- 
pact way to express these LHV violations. Replacing 
Pupb with p SC p, all 4 inequalities (@J change sign and the 
contradictions disappear. The argument used is of the 
same type of |lfj). Since we have mixed states, we cannot 
formulate it rigorously in terms of "signs of the eigenval- 
ues" of wavefunctions that are eigenfunctions as in the 
formulation of 0,0]. However, assumption (5) of is 
here verified and guarantees that each of the observables 
we consider (i.e., all those with nonzero expectation) can 
be traced back uniquely to one of the 4 ket states 

Preparation. The second question we are interested in 
is how to create the state pupb- While p sep is a mixture 
of separable pure states and can be produced by standard 
methods, the reflection operation Se4 is nonunitary and 
no quantum circuit is known for it. It is obvious that if 



7] 
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we start with any separable state and consider the orbit 
given infinitesimally by the Lie algebra Lie{— iad\ jkl }, 
j, k, I, — 0, 1, 2, 3, this will certainly contain all states (in- 
cluding entangled) and only states. In particular, it is 
possible to find explicitly the Hamiltonian of a unitary 
transformation that maps p scp to pupb- This is the fol- 
lowing concatenation of nonlocal constant infinitesimal 
generators: 



-zadA 3 



for t e [0, 



— iadn = 



-iad 



(A011+A033+A101 +A110+A303+A330) 



. (5) 

where t p = 2v / 2tt is the period in the parametriza- 
tion we are using. The effect of the first part is to 
switch the sign of the entire homogeneous tensor £, ex- 
cept for the 6 2-coherences with pairwise equal indexes 
e {on,033,ioi,no,303 l 330} < In spite of the non i OC ality of 

cxp (— itadA 333 ), its action on p sep has remarkable "local- 
like" properties at the end of the time interval t = 5- . In 
fact, it changes sign to all and only the components of 
index "1" in each of the 3 qubits of the 4 components of 
p scp . This can be understood by means of the tensored 
Rodrigues' formula presented in |12j . which provides an 
explicit closed form expression for one-parameter flows 
such as each piece of (JSJ. For example: 



cxp HiadA 3 3 3 ) = Efclo L -^~ 



adA 



r®3 



£\/2sin(^=)ad A 3 33 - 2 ( 1 - cos(J=) ) ad 



V2- 



] 2 

A333 I 



where ad A = ^(ad^, ® aad^ 



aadx„+ aadl <g> adx 



aadi 



aad> 



aad> 



ad? 



ad 



ad 



ad 



with ad A , 3 = V2i(5 32 - 823), ad^ 3 = 2(S 2 2 + 833), 
&&d\ 3 = \/2(<5i 4 + &11), aad^ 3 = 2{5\\ + 844). Call- 
ing p^t = exp (-j^adA 333 ) p sop , then we have p int = 
lEi=il^)<^|with|^) = |01-),|l-0),|-01),|+++). 
Hence p m t is still a separable state with the same prop- 
erties as Psep- The second piece of (0 is the one crucial 
for the purposes of creating the entanglement. Its action 
can be analyzed by means of arguments similar to the 
ones above. Both pieces of (J2J) are nonlocal and they 
must be so because of the randomness of all 1-qubit 
densities. While at the begin and end of each interval 
the density is PPT, during both evolutions there is 
always bipartite entanglement through each 1-2 cut. 
Exchanging the order of application of the two Hamilto- 
nians of J5J, one still gets a map p sep — > pupb- However 
the intermediate state reached is different: p lnt = 

e_xp ( — «X a d(Aoii+Ao33+Aioi+Ano+A303+A33o)) Psep = 

S , 64(/°int)- It is straightforward to check that also pj nt 
obeys to (QJ, hence it is bound entangled. 

A PPT entangled orbit. Perturb pupb by applying a 
Hamiltonian for a certain time interval. From random- 
ness of all the 1-qubit reduced densities, we have that 



Pupb is (trivially) invariant to all local actions, as al- 
ready mentioned above. When a 2- or 3- spin Hamilto- 
nian is applied, the resulting density is generically not 
PPT. There are at least two exceptions: the directions 

— zad(A 011 +Ao22+A 33+Aioi+Aiio+A 202 +A 22 o+A 30 3+A 330 ) and 
— iadA 222 . While the first generator induces no action 
at all on pupb (he., pupb is a fixed point for it), along 
the one-parameter orbit of the second one the density 
matrix is always PPT with respect to any bipartite cut. 
Each density operator on such orbits has rank 4. This 
orbit does not modifies any of the 6 1- or 2-qubit reduced 
densities, but affects only the 8 3-qubit coherences: Q ln , 
g 113 , g 131 , g 133 , g 311 , g 313 , g 331 , g 333 , i.e., exactly those 
3-coherences that certainly do not yield bipartite entan- 
glement (no index "2", see 8]). We can look at what 
happens for p or b(i) = exp (— iiadA 222 ) p S ep and then con- 
sider the corresponding reflection p or b(i) = >§64 (Porb(i)) 
[l4|. Computing the sum of the series as in ©, 

Porb(i) = Psep - Sin(-^=)a;A { ii3 ) i3i ) 3i li 333} 



1 - C0S ("^f)) ^ A {111, 133,313,331} (7a) 



Q {a} ^{a} - Sm(-^=)xA {11 3 43l! 3 11:3 33 } 
COs(--^)a;A{ 11 l,133,313,331} 5 



t 



(7b) 



(6) where {a} denotes the subset of indexes jkl con- 
taining all 0- and 2-coherences of p S e P : a = 
000, 011, 013, 031, 033, 101, 103, 110, 130, 301, 303, 310, 330. 
In words, on the one-parameter orbit J7J) all the 

2- coherences are constants of the motion, while the 3- 
coherences evolve according to a sinusoidal law. Neither 
of the last two terms alone in (|7b|) is a density (they are 
both traceless). The sinusoidal law in JJJ has the effect 
of "swapping" the bound entanglement between p or b 
and p orb . In fact, if p O rb(0) = p scp and p orb (0) = pupb, 
at t = -J, sm (jT/|) = 1 and we have to replace the 4 

3- coherences of Q and with the other 4. Denoting 
Poq = Porb (if): 



Poq 



,000 _ 1 

2\/2 



{011,013,101,110,130,301} _ x 
{031,033,103,113,131,303,310,311,330,333} _ 



cP kl = otherwise. 



(8) 

It is easy to check that in this case p rb(-f ) i s separable 
and given by porb(x) = 3 £*=i Poq 3 : 
Po qi = (g°X + g 1 X 1 ) <g> (g°\o ~ g 3 X 3 ) <g> (g°X - g 1 Xi) , 
Poq 2 = {g°X - 0%) <8> (g°X + g 1 Xi) <g> (g°X Q - g 3 X 3 ) , 
Poq 3 = (g°X - g 3 X 3 ) ® (g°X - 0%) (g°X + Q 1 Xi) , 
Poq 4 = (g°X + g 3 X 3 ) <g> (g°X + g 3 X 3 ) <g> (g°X + g 3 X 3 ) , 

(9) 
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or, in terms of kets, p oq 



j)(6j\ with \0j 



•1- 



From Q, the conjecture is that both p OI b{t) and p rb(t) 



+1), |1 - +), 1 000), which is another UPB basis. Hence 
p oq (obtainable from J2J) is bound entangled. The explic- 
itly known mixture of the separable state p or b(if ) tells 
us that the triplets of observables to be used to detect 
violations of LHV similar to Q in p oq are now different: 



tr (/Oo q A 3i) tr (poqAioi) tr (p O qA 130 ) < 0. 
tr (poqAois) tr (p oq Ai 3) tr (p qAno) < 0. 
tr (poqAon) tr (p qA 3 oi) tr (poqA 3 i ) < 0. 
tr (p oq Ao33) tr (poqA 3 03) tr (poqA 3 3 ) < 0. 



(10) 



The explicit knowledge of the mixture of pure states is 
crucial in setting the inequalities (@J and (|10|l . Notice in 
fact that Pupb is compatible with the LHV model given 
by the triplets of (fTU|) and, similarly, that p or b(x) * s 
compatible with those of Q. 

After another fourth of period, the orbit p or b(i) reaches 

another separable state p O rb(-f0 = J Ej=i \4>j)(4>j\: with 
\<j>i) = |10-),|0-1),|-10),| + ++), or, in terms ofBloch 

^4 



vectors, porb(if) = \ El 



. Poll. 



with 



A*! = 
Poh 2 = 
Poh 3 = 
Pohi = 



(g°\ - g 3 X 3 ) 
{g°X + g 3 X 3 ) 
{g°X - g'Xi) 
[q°Xo + q 1 Xi) 



(g°X Q + g 3 X 3 ) 
(g°Xo - g^i) 
(g°\o - g 3 X 3 ) 
(g a Xo + 0%) 



(g°X - ^Ai) , 
(g°X - g 3 X 3 ) , 
(g°X + g 3 X 3 ) , 
(g a X + e x Ai) . 



Since cos(^^) = — 1, Porb(^f-) has the same tensor 
as Pupb, but with the opposite sign in the 3-coherences. 
Comparing p sop and p or b(^) (see Fig. QJ, p or b(if) cor- 
responds to a change of the original UPB basis in which 
each member of the ket \<pj) is obtained by means of a 7r 
rotation on the corresponding Bloch sphere around the 
A2 axis. For example, \ipi) = |0) <g> |1) ® |+) rotated 
around A222 = A2 <8 A2 <8> A 2 gives \(f>i) = |1) ® |0) ® |— ) 
(first row in Fig. and so on. A similar relation holds 
between p or b(x) an d Porb(-p-)- The effect of the nonlo- 
cality of the rotation applied is the swapping of bound 
entanglement between p or b and p or b, plus the recombi- 
nation of pieces of product states into different (rotated) 
product states. While the "blocking" state p SC p 4 under- 
goes a rotation around A2 on each qubit (plus "swap- 
ping") to yield p oq4 at t = ^ and p oh4 at t = ^ (see 
last row of Fig. [IJ, every other triplet of Bloch vectors 
of Poq comes from a ? rotation of 3 Bloch vectors of 
p scp belonging to 3 different pieces Pscp^, j — 1,2,3. For 
example, p oqi = | + 1— )(+l — | (i.e., the squares in the 
first row of Fig. ^) contains the ^ rotations of the first 
Bloch vector of p SC p 1 , the second of p S op 2 and the third 
of Psep, (he., the f rotations of the diagonal bullets in 



Fig.HJ. While at t 



state admits a neat splitting into separable and bound en- 
tangled, the situation is more ambiguous when t 7^ -r 2 -. 




FIG. 1: The Bloch vectors of the 4 pure states of the separable 
densities p scp = p or b(0) (•), p rb(x)( c ) and porbC^f) (A). 



share a percentage of separability and entanglement in a 
"nonseparable" way. There may be other nontrivial PPT 
entangled orbits, but they are more difficult to find. 

Finally notice that, as a byproduct, one gets a 
simpler scheme for the preparation of pupb from a 
separable state: pupb = exp(rfadA 222 )p or b("J) = 

3r„ 



exp ( - ited A222 ) p orb ( -f ) 
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